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Abstract. 

The purpose of this paper is to study the thermodynamic equilibrium properties 
of a collection of non-interacting three-dimensional (3D) magnetically anisotropic 
nanoparticles in the light of classical statistical physics. Pertaining to the angular 
dependence (a) of the magnetic field with the anisotropy axis, energy landscape plots 
are obtained which reveal a continuous transition from a double well to a single well for 
a = ^ and show asymmetric bistable shape for other values of a. The present analysis 
is related with the interpretation of equilibrium magnetization and static susceptibility 
of nanomagnetic system as a function of external magnetic field, B, and temperature, 
T. The magnetization and susceptibility confirms the non Langevin behaviour of 
magneto-anisotropic monodomain particles. The susceptibility analysis establishes the 
ferromagnetic, antiferromagnetic and paramagnetic like coupling for various a. This 
study reveals the essential role of magneto anisotropic energy in the interpretation of 
the magnetic behaviour of a collection of noninteracting single domain nanoparticles. 

PACS numbers: 75.75.H-a, 75.50.Tt, 75.40. Cx,75.20.-g 
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1. Introduction 

Nanometer sized magnetic particles have provoked immense interest in both scientific 
as well as technological arena [H [21 [3], [32], [5]. The development of intense fabrication 
techniques helps the preparation of nanoparticles with satisfactory structural and 
chemical properties. The study and analysis process gains acceleration due to the 
exaggerating growth of measurement facilities like magnetic force microscopy [6] , micro- 
SQUIDS [7j and other magnetometry measurements P,|9]. Such techniques have lead to 
the measurement of the magnetization process of single magnetic clusters in nanometer 
scales. 

The magnetic moment of the nanoparticle consists of single domain structure of 
ferromagnetic spins with a large net spin, S (~ 10^ — 10^) and hence it is named as 
supermoment [IDl [HI [12] • This spin couples with a large number of environmental 
degrees of freedom of the host material. Dynamical disturbances of the surrounding 
environment leads to a rotational Brownian motion of the large spin surmounting the 
magnetic anisotropy potential barriers [131 [T^. the high barrier limit, the magnetic 
response of the non-interacting single-domain particles follow the Neel relaxation process 
with the relaxation time r characterized by the relation 



where Tq ~ 10^^'^ — 10^^"^ sec. Here Tq is related to the intra- well motion and the height of 
the energy barrier due to anisotropy is AEa = KV where K is the anisotropy constant, 
V is the particle volume, fc^ is the Boltzmann constant and T denotes the absolute 
temperature. Depending on the relation of r with the measurement time tm-, various 
interesting phenomena can be observed. For t tm, the magnetic moment exhibits 
the thermal equilibrium distribution of a paramagnet. For r ^ tm, the magnetic mo- 
ment stays very close to the energy minima as the reversal mechanism is blocked. For 
r ~ tm, nonequilibrium phenomena i.e. magnetic relaxation is observed. In this work, 
all discussions are concentrated in the thermal equilibrium regime where r <^ t^. In 
this context, the classical Stoner-Wohlfarth (SW) simplistic uniform rotation model of 
giant magnetic vector really provides insightful and realistic picture of magnetization 
reversal of single domain nanoparticle ^l5j. Here, we follow the same kind of arguments 
as proposed by Stoner and Wohlfarth [T5] . 

The first study on the superparamagnetic behaviour of an aligned assembly of 
uniaxially anisotropic particles was made by F. G. West [16]. Further investigations 
regarding possible configurations encountered in experiment was made by Miiller and 
Thurley[T7]. Deviations from classical Langevin theory was demonstrated by several au- 
thors [m [ini [201 [211 [221 [23] . Madsen et al have studied the effect of anisotropic energy on 
the interpretation of magnetization data for antiferromagnetic particles [24J. Vargas et al 
depicts a second order phase transition in non-interacting magneto-anisotropic nanopar- 
ticles when the external magnetic field is applied perpendicular to the anisotropy axis 
with the order parameter being the magnetization parallel to the field [25l [26] . Magneto- 
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caloric properties of noninteracting magneto-anisotropic nanoparticles have been studied 
in [30]. 

In this work, we investigate the effect of magneto-anisotropic energy on the equilib- 
rium thermodynamic properties of fine particles. By thermal equilibrium behaviour, we 
meant that the measurement or observation time, tm, is much larger than the character- 
istic relaxation time, r, of the system i.e. we restrict our discussion in the regime t^n ^ t. 
On the other hand, below a certain critical size (for Fe, it is 150 nm.), it is not energet- 
ically favourable to form a domain wall, and the particle is said to be monodomain or 
fine particle system [11]. Having explained the terms "equilibrium properties" and "fine 
particles", we are now ready to describe briefly about our findings. Variations in the 
angle between the anisotropy axis and the external magnetic field, a, and between the 
former and the magnetic moment, 6, give deep insights into the thermodynamic equilib- 
rium properties of magneto-anisotropic nanoparticles. We have extended the study of 
Vergas et al [251 126] by choosing the magnetic moment and the anisotropy field vectors 
to be independent and arbitrary. On the other hand, Vergas et al, restricted their study 
for the particular case in which external field is perpendicular to the anisotropy axis. 
Also, we consider the particle size distribution of nanoparticles rather than the case of 
identical noninteracting particles as studied by Vergas et al [251 126] . Thus the present 
study is much more realistic and close to the experimental realizations [271 [281 129]. 

The variation of the magnetization and susceptibility with respect to external pa- 
rameters like external field, 5, and temperature, T, is extensively studied. The effect 
of anisotropy is evident from the magnetization versus reduced magnetic field, ^, curve. 
The variation of inverse susceptibility with temperature for different angles, a, shows 
paramagnetic, ferromagnetic and antiferromagnetic like coupling which clearly exhibits 
the effect of mgneto-anisotropic energy on the static susceptibility of nanomagnetic sys- 
tem. 

With the preceding background, the rest of the paper is organized as follows. In 
the next section, we discuss about the model system and other basic considerations 
about this model system. Section 3 deals with the model Hamiltonian and the relevent 
thermodynamic functions. In subsection 4.1, the energy landscape for the nanomagnetic 
system as a function of a is explored. The variation of equilibrium angle, 9, with di- 
mensionless reduced magnetic field, h, for different a is demonstrated in the subsection 
4.2. Section 5 deals with the variation of magnetization and susceptibility with respect 
to ^ and T. Finally, the paper is concluded in section 6. 

2. Model and basic considerations 

In this section, we discuss about our model and some basic considerations which one 
needs to study thermodynamic equilibrium properties of single domain particles. Kittle 
has shown that below a critical size, domain wall formation is energetically unfavourable 
[TT] . This kind of particles are called single domain particles. In the absence of an ex- 
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ternal magnetic field, a bulk ferromagnet may have no net magnetization due to the 
cancellation from different domains. But, a single domain particle acts as a giant mag- 
netic moment and the magnetic moment per particle depends on the particle volume 
and the number of atoms it has. Throughout this paper, we concentrate on the thermo- 
dynamic properties of a collection of such monodomain particles which are dispersed in 
solid matrix [311 132] ■ Further, we restrict our study to mathematically tractable systems 
with axially symmetric magneto- anisotropy. But, it provides valuable insights into more 
complex situations. All our considerations are based on the Stoner-Wohlfarth model for 
single domain particles [15]. In this model, all spins within the particle are aligned due 
to exchange interaction and this is the dominating magnetostatic effect within the par- 
ticle. The giant magnetic moment direction fluctuates, because the anisotropy energy 
is comparable to thermal energy. The direction of the moment is determined by the net 
anisotropy of the system and energy is minimized. This magnetization reversal process 
occurs by coherent rotation i.e. the atomic spins remain parallel to each other as they 
rotate to a new direction [T3l [H]. Here, we are considering ideal monodomain parti- 
cles in which other more complex interactions both within the particles and between 
the particles are neglected. Thus, our system has only magnetocrystaline energy and 
Zeeman energy due to the interaction with external field. 

Now, we need to specify the validity of this model, especially in which temperature 
range it is valid. To specify this temperature range, we follow the arguments given 
by J. L. Garcia Palacios [20]. It has already been mentioned that the thermal equi- 
librium behaviour of ideal monodomain particles with uniaxial anisotropy is observed 
when tm ^ T. From equation (1), one can easily understand that thermal equilib- 
rium behaviour exists when anisotropic potential barrier, AEa, much larger than the 
thermal energy, ksT. Besides, the "high-barrier" regime, equation (1) still holds down 
to AEa/ksT > 2. It is known that tq for magnetic nanoparticles is ~ 10^^'^ — 10^^^ 
sec. Thus, the thermal equilibrium range for a given measurement time, tm, is given 
by In > AEa/ksT > 0. For magnetic measurements, t^ ^ 1 — 100 sec and the 
thermal equilibrium range is quite wide 25 > AEa/kBT > 0. So, the frequently en- 
countered statement that the thermal equilibrium behaviour occurs when AE^ > ksT 
is needlessly constrictive. For example, if tm = rol0^2 sec (typical value for magnetic 
measurements), one finds that AEa/ksT ~ 27.6. For AEa/ksT = 25, one can obtain 
r = O.OStm- Thus, the system is in thermal equilibrium, but AEa is still much larger 
than ksT 

3. Hamiltonian and other relevent quantitites 

We consider a collection of non-interacting single domain magnetic particles. In a single 
domain particle, nearly 10^ magnetic moments are coherently locked together in a given 
direction thus yielding a supermoment. In this nanomagnetic system, every particle 
consists of a single magnetic domain with all its atomic moments are rotating coherently 
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and resulting in a constant absolute value of magnetization m = mgV, where V is the 
volume of the particle and is the saturation magnetization which is supposed to 
be independent of particle volume and temperature. In our system, the Hamiltonian 
consists of two parts, one representing the Zeeman energy and the other is the anisotropic 
energy (due to the crystalline structure of the particle). For the sake of simplicity, we 
consider a temperature independent uniaxial anisotropy. Denoting the external applied 
magnetic field as B, the Hamiltonian is given by 

Tiim) = — (m-n) —rn-B, (2) 

V / 

where n is a unit vector along the anisotropy axis, m is the magnetic moment of the 
single domain particle and B is the direction of the external magnetic field as shown 
in Fig. 1 and K is the anisotropy constant. This model is valid only if the exchange 




the angles a, 9, and (p as referred in the text. 



interaction strength of the system is much larger than K and B. Now denoting {6, 0) 
and {a, 0) as the angular co-ordinates of rn and B respectively and choosing n as the 
polar axis of the spherical polar co-ordinate system, one can write the total magnetic 
potential in the form as follows, 

— PH = a cos^ ^ + ^11 cos 9 + ^± sin 9 cos (f), (3) 

where a = = ^cosa, = ^sina, ^ = (3 = ks is the Boltzman 

constant and T is the temperature of the system. One can rewrite Eq. (3) as follows 

Ti-eff = = cos^ 9 + 2/i(cos a cos 9 + sin a sin 9 cos 0), (4) 

a 

with h = Expressions for the magnetization curve where anisotropy is included have 
been discussed by [iHl UHl EQl [21]. This calculation is based on the theory of Hanson et 
al [18] and Respaud [21] . For uniaxial anisotropic nanomagnetic system, the anisotropy 
energy depends only on the angle between the magnetization and the direction of easy 
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magnetization in the particle, 9. In thermal equilibrium, the magnetization in the 
direction of in is proportional to the Boltzmann factor for a fixed orientation of the easy 
axis : 

f{rn) — z~^ex^ a{m ■ + ^{rh ■ h) , (5) 

where m, h are the unit vectors along the direction of magnetic moment and the external 
magnetic field and z is the partition function defined by 

{m-hf + 2h{m-h)]. (6) 



z = I exp 
Thus, we have 

m-h — cos A = sin q; sin 9 cos + cos a cos 9 (7) 
m-n — cos 9. (8) 

Using the Boltzmann statistics, the expectation value of the reduced magnetization with 
a given orientation of the easy axis is given by 

M /o'" #/o"cosAe-^«// sm9d9 

mia) — — = (cos A) = — — Ty- — - — - — — . 9 

^ ' rus ^ ' J^''d(f)J^e-^^ff sin 9d9 ^ ^ 

None of the integrations in Eq. (9) is doable analytically. However, equation (9) can be 
simplified by using the modified Bessel functions and performing the analytic integration 
over (f). Thus, the magnetization of such a collection of non-interacting identical particles 
aligned with an angle a with respect to B is given by 

m(a) = ^. (10) 

with 

N(a) = / d9 sin 9 exp(2h cos A + cos^ 9) [sin a sin 9Ii(2h sin ct sin ^) 
Jo 

+ cosQ;cos^/o(2/isinQ;sin^)], (11) 

/"IT 

D{a) = / (i6'sin6'exp(2/icosA + cos^6')/o(2/isinasin6'), (12) 

JO 

where Jq and Ji are the modified Bessel function of order and 1 respectively. 
Considering random distribution of anisotropy axes one can show 

Msih) — - dasinam{a). (13) 



2 

However, there will be a distribution of particle sizes in any real fine particle system. 
The existance of particle size distribution can be taken into account by taking average 
over the full particle size distribution. Thus, the magnetization of such a system 
with a distribution of particle sizes is consists of the sum of contributions from the 

superparamagnetic and the blocked particles. The weightage of these two is maintained 
by the size distribution function of the particles, f{y)- Now, the magnetization for this 
polydispersive system is given by 

rVsp /-oo 
Mpoi= M,p(y,h)f(y)dy+ MB(h)f(y)dy, (14) 

Jo •'Hap 
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where y is the reduced volume ^ with the mean volume Vq and ysp = Vgp is 

the critical volume for superparamagnetism which is given by = ^^j^- Msp{y,h) 
and Msih) are the reduced magnetization for the superparamagnetic and blocked 
particles respectively. It is known that Msp{h, y) = msL{2ayh) where Langevin function 
L{x) = coth(a;) — ^. Usually such anisotropic nanomagnetic system follows log-normal 
distribution of particle size i.e. 

/(y) = ^exp[-<i|^], (15) 



where 7 is the dispersion of the corresponding distribution. For the numerical 
integration, we have used the following parameters : 7 = 0.8, Vq = exp(— 7^/2), and 
average blocking temperature, (Tb) = 15.5K, [33]. T is measured in units of (Tb), and 
h is measured in units of m^-B / K. 

Numerical integration programme in FORTRAN are performed to calculate 
magnetization for the monodispersive and the polydispersive system by using equation 
(13 ) and (14) respectively. Static magnetic susceptibility of the polydispersive system 
is defined as 

Xpoi = -g^- (16) 

Now, we have defined our system and other essential thermodynamic functions. In the 
next two sections, we analyze the thermodynamic behaviour of such magneto-anisotropic 
nanomagnetic system. 



4. Energy barrier and equilibrium angle 

In this section, we discuss about the behaviour of the magnetic potential energy as a 
function of the several parameters used in the Hamiltonian of the present nanomagnetic 
system. The variation of equilibrium angle, 6, between magnetic moment and easy axis 
of magnetization as a function of h is also studied in this section. 

4.I. Energy barrier 

One can rewrite equation (4) as follows 

U {6, 0) = — — ^ ' = sin^ 6 — 2h{cosacos6 + sin a sin 6 cos 0), (17) 



a 



where (3 = -j^. The stationary points for Equation (17) occurs for = and 



TC. 



The stationary point for = vr corresponds to a maximum, so it is of no physical 
interest. On the other hand, stationary point = corresponds to a maxima at 6m and 
minima at 61 and 62- One can determine two equilibrium directions of the magnetization 
associated with polar angles 61 and 62 (lie in the x-z plane) from the following condition 

f^0,^>0. ,18, 
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Figure 2. (color online) Energy is plotted against the angle 6 for different cases with 
: (a) a = = 0; ft. = 0.25 (black circle), and h — 1.1 (blue square) (b) a ~ j (p = 0; 
h = 0.1 (black circle), and h = 0.6 (blue square) (c) a = ^ (j) = 0; h — 0.7 (black 
circle), h — 1.01 (blue square) {d) a — tt tj) — 0; h — 0.25 (black circle), and h — 1.1 
(blue circle). 



and the saddle point is determined by 
dU 



0,^^ < 0. 



(19) 



On the other hand, one can determine the critical value of the ratio of field to barrier 
height {he) at which the potential loses its bistable character by using the following 
condition 



de 



Using equation (20), one can easily show that 

1 



2 2 

COS3 a + sins a 



One can rewrite equation (21) as follows 

27, 



h 



2\3 



-h'^. sin^ 2a 
4 



0. 



(20) 



(21) 



(22) 



It can be shown that lies in the range 0.5 < \h(.\ < 1. h^. = 1 occurs for a = 
or a = |. Whereas h^ = \ can be seen for a = j. It is very unlikely that one 
can derive all the derivatives and hence the barrier heights {Bi and B2) for arbitrary 
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a. It is easy to derive barrier heights and to know the nature of the potential for 
some particular values of a. For a = 0, -Bi = (t(1 + and B2 = cr(l — hy. Thus 
the potential has the asymmetric bistable form as shown in Fig. 2(a). On the other 
hand, this asymmetric potential becomes symmetric and the barrier height becomes 
Bi = B2 = a{l — hY for a = |, = (see figure 2(c)). Again for a = f , the potential 
becomes asymmetric bistable below h = 0.5 (see figure 2(b)) and the barrier heights be- 
come [22j El = 2(7^(1/2) - (/iV2) - h^'^/A) + {l/2){^ {K^ / A) + (1/2) - 3/i/2) and 

B2 = Bi/2 + a^{l/2) - (/i72) + h^i^K'/^) + (l/2)(^(/iV4) + (1/2) + 8/1/2). For 
a = TT, one can show that Bi = a{l — hY and B2 = cr(l + h)"^ and is just the opposite 
of the case for a = (see figure 2(d)). In general the potential retains its asymmetric 
bistable character for < h < and a ^ ^. 

It is evident from fig. (2) that the potential energy of this nanomagnetic system has 
two minima separated by a maxima if and only if h is less than a ceratin critical value 
(he) which varies from 0.5 for a = | to 1.0 for a = 0, |. So, the system consists of two 
potential barriers Bi and B2 which are in general unequal except for the case of a = ^. 
li h > he the bistable character of the potential disappears and the system has only a 
single maxima or a single minima. At h = h^., the second minima becomes a point of 
inflexion which is clearly seen in figure (2). 



4-2. Equilibrium Angle 

Here, we are considering a collection of non-interacting magnetic single-domain nanopar- 
ticles in the presence of an external magnetic field. The variation of the equilibrium 
angle between the anisotropy axis and the magnetic moment, 6, versus h for such a 
collection of noninteracting nanomagnetic system is shown in Fig. 3. Now, using the 
maxima and minima condition of U{6, (p) as defined by equation (18) and equation (19), 
one can obtain the equilibrium angle for the magnetization direction in the x-z plane for 
some specific values of a. Finally, the expression for the equilibrium angle 6 for = 
under different orientations of a are as follows : 

(i) for a = and = 0: 

= cos-^(-/i),0i,2 = 0,7r (23) 

i.e U{6, 0) has minima at 61^2 = 0, tt and maxima at 6m = cos^^{—h). 

(ii) For a = f and = 0, 

the minima of the potential are at ^ = ^1 = ^ - sin"^ ( ~ | + ^^2^^ ); < ^1 < and 
at e = 02 = ^ - sin-^ (I + ^^^^); X < ^2 < TT. On the other hand, the magnetic 

potential has maxima at 6 = Om = j + sin"^ (| + ^^^); f < < x 133!- 

(iii) For a = | and = 0, 

U{9,(f)) has minima at = ^1 = sin~^(/i) and 9 = 62 = n — sin^^(/i) and maxima at 
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(iv) and finally for a = n and (f) = one can show that minima of the magnetic potential 
are at = 612 = 0, vr and maxima at = = cos~^{h). 

In Fig. 3, we plot these solutions as a function of h. Figure 3(a) shows that 6^ for 
a = 0; = increases from to 3.0 as h is altered from —1 to +1. On the other hand, 
9m for a = 7i; (p = decreases from 3.0 to as is varied from —1 to +1 (see fig. 3(b)). 
From figure 3(c), it is evident that as h is varied from —1 to +1, 61 increases from — f to 
I (black filled circle) and 62 decreases from ^ to | (blue filled square). It is seen from 
fig. 3(d) that 9i (black filled circle) and 6^ (blue filled square) increases monotonically 
to a maximum value of ^ and ^ respectively at h = 0.5; on the other hand, 6*2 decreases 
monotonically to a minimum value of ^ at /i = 0.5 (red filled triangle). 



5. Magnetization & Susceptibility 

Magnetization and susceptibility are the most fundamental thermodynamical quantities 
of non-interacting magnetic nanoparticles with axially symmetric magnetic anisotropy. 
In this section, we analyze about the variation of these two fundamental quantities with 
temperature and externally applied magnetic field. The differences and similarities of 
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the magnetization and susceptibility between the ideal superparamagnetic system and 
a colletion of noninteracting anisotropic monodomain particles are presented here. 

5.1. Magnetization 

The magnetization along the direction of the external magnetic field for classical spins 
with axially symmetric magnetic anisotropy is defined by equation (13). We illustrate 
this magnetization as a function of ^ using equation (13) for a system of identical non- 
interacting monodomain particles for different values of a in Fig. 4. One can easily 



1 



0.8 
0.6 

CD 

0.4 
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2 4 6 8 10 

Figure 4. Reduced magnetization for the monodispersive system as a function of ^, 
for various values of the anisotropy parameter, a. The black circle and blue square 
represents Langevin (cr = 0) and Ising {a » 1) cases respectively. Green square 
and red upward triangle represents magnetization curves for a = 5.0 and a = 2.0 
respectively. 

observe that magnetization curves differ from the Langevin law in all the cases. As a 
decreases, the difference vanishes and with the increase of a, the magnetization curves 
become closer to the Ising case. In the limits of high and low field, ^, magnetization, 
Mb, approaches the Langevin value. The largest influence of anisotropy can be ob- 
served in the intermediate field regime. This confirms the non-Langevin behaviour of 
a collection of noninteracting magneto-anisotropic single domain particles. In Fig. 5, 
we show the variation of magnetization with reduced temperature, = jf^, for two 
different values of h. Here (Tb) is the average blocking temperature of the system |35j . 
For convenience, we also plot the magnetization curves for the Ising and Langevin cases. 
It is clear that magnetization versus temperature curves for the anisotropic magnetic 
nanoparticle system show maximum near = 1.0 unlike the Ising and Langevin cases. 
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Figure 5. Reduced magnetization along the external field axis for a collection of 
noninteracting, monodomain and polydispersive nanoparticles as a function of reduced 
temperature. The black circle, blue square and pink upward triangle represents 
Langevin, Ising and magneto-anisotropic cases respectively: (a) h=0.2 and (b) h=0.5. 

This maximum can be interpreted as follows. From equation (14), one can observe 
that Mpoi has two parts, the superparamagnetic contribution and the blocked particle 
contribution. As the temperature increases the fraction of the superparamagnetic con- 
tribution increases till temperature reaches the blocking temperature. Now, above this 
blocking temperature at which maximum in magnetization is observed, the system usu- 
ally becomes superparamagnetic and magnetization decreases rapidly with the increase 
of temperature due to the thermal agitation. Thus, in equilibrium and for h < he, 
one can observe maximum in magnetization for the polydisperse magneto-anisotropic 
nanoparticle system. One can observe that the maximum of the peak is exactly not 
at Tj. = 1.0, but with the increase of field it shifts to lower relative temperature i.e. 
maximum are seen at = 0.89 and at Tr = 0.83 for h = 0.2 and h = 0.5 respectively. 
Also one can notice that as the field increases, the peak becomes broader. 

5.2. Susceptibility 

In order to demonstrate the effect of magnetocrystalline anisotropy on the thermody- 
namical quantities, we demonstrate the behavior of susceptibility of the non-interacting 
nanoparticles in presence of an external magnetic field. We plot the susceptibility and 
inverse susceptibility curves with respect to temperature for different orientations in a 
between the anisotropy axis and external field in Fig. 6 and Fig. 7 respectively. Again, 
susceptibility shows a maximum at finite temperatures (T^ = 0.89 and = 0.83) 
for h = 0.2 and h = 0.5 respectively. Also, as the field increases the peak of the 
susceptibility curves become more broader. In order to demonstrate the effect of mag- 
netocrystalline anisotropy, we illustrate the behavior of inverse susceptibility of the 
noninteracting nanoparticles with reduced temperature in the presence of an external 
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Figure 6. Static magnetic susceptibility along the external field axis for a collection of 
noninteracting, monodomain and polydispersive nanoparticles as a function of reduced 
temperature. The black circle, blue square and pink upward triangle represents 
Langevin, Ising and magneto-anisotropic cases respectively: (a) h=0.2 and (b) h=0.5. 

magnetic field for different orientations of a between tlie anisotropy axis and external 
field in figure 7. For a = 0, = (in black circle), one observes a susceptibility curve 
resembling a system with ferromagnetic like interaction and it follows Curie- Weiss law. 
The same kind of ferromagnetic like coupling is prevailed for a = (in orange square) 
and for a = ^ {in pink upward triangle). This means the zero crossing at tempera- 
ture axis occurs at positive values. When the angle between the anisotropy axis and 
magnetic field is a = ^ (in red downward triangle), the system resembles paramagnetic 
behavior and it follows simple Curie law. However, for a = | (in blue diamond) and 
for a = I (green filled square), a Curie- Weiss antiferromagnetic-like behavior is ob- 
served. By extrapolating one can observe the zero crossing in the temperature axis at 
negative values. Thus, the extrapolating transition temperature continuously changes 
from positive values to negative values. From this kind of observation, one can conclude 
that the anisotropy field acts as a ferromagnetic, antiferromagnetic or paramagnetic like 
coupling among the magnetic nanoparticles depending on the relative angle between the 
anisotropy axis and the external field. The net effect of the anisotropy field resembles a 
continupus transition from a ferromagnetic to an antiferromagnetic like coupling. 

6. Conclusions 

In conclusion, we have confirmed that a collection of noninteracting magneto-anisotropic 
particles cannot be described with the help of the classical Langevin theory, i.e., their 
thermodynamic equilibrium magnetization and static magnetic susceptibility cannot be 
described by the Langevin function and its derivative. This deviation is due to the 
presence of magneto-crystalline anisotropy. The effect of magnetocrystalline anisotropy 
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Figure 7. Inverse magnetic susceptibilities as a function of reduced temperature. 
Different curves represent different angles between anisotropy axis and external field 
(for details see text). 

is explored through magnetization curves and susceptibihty curves. The variation of 
inverse susceptibihty with temperature shows paramagnetic, ferromagnetic and anti- 
ferromagnetic couphng behaviour for different orientations of a. We also present a me- 
chanical analogy for a system in 3D as a frictionless particle moving on a sphere rotating 
about its vertical diameter. This study reveals the essential role of the magnetocrys- 
talline anisotropy energy in interpreting equilibrium magnetization and susceptibility of 
a collection of non-interacting single domain nanomagnetic particles. 

Appendix 

In this appendix, we discuss about the mechanical isomorph of the nanomagnetic system. 
We consider a rigid sphere of radius Rs rotating along its vertical diameter at angular 
frequency u with a frictionless particle of mass M free to move on the surface of 
the sphere. In the rotating co-ordinate frame (r, 6, (p) attached to the rigid sphere, 
in addition to the gravitational force on the mass Mcfa = MgaCosOir — Mga sinOee, a 
fictitious centrifugal force / is acting on the particle and is given 

Mu'^r sin 6 cos 6eg. An external field F is acting on the azimuthal plane such that 
F = Fsin^e^. Now, the kinetic energy of the particle is given by 

T = ^M(r'^ + r^e'^ + r'^sm'^e^^y (24) 

The potential energy of the particle consists of three parts U = Ug^ + Uc + Ue, where 
Ug^ is coming from the gravitational force field part, Uc is the fictitious centrifugal part 
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Figure 8. Mechanical rotating system analogous to nanomagnetic system 



and C/g is the external field part. Now, one can easily find out the three parts of the 
effective potential energy as follows : 

Ug^ = -Mg^R, cose 

Ue = FRsSme COS (p (25) 
Thus the Lagrangian of the system becomes 

C = ^M(f^ + r^^^ + sin^ ^0^) + MgaRs cos 9 
1 1 

+ -Muj'^R] - -Mu^ cos^ 9 - FRs sin 9 cos 0. (26) 
Now, the effective potential energy which includes the effect of gravity, rotation of the 

— * 

system and the external field force F is given by 

U , UJ^Rs u'^Rs Or. F ■ . 

u — — = — cost^ 1 cos 9^ smt^cosffl. (27) 

MgaRs ^ga Mga 

Comparing equations (3) and (27) one can easily understand the analogy between the 

mechanical system and the magnetic nanoparticle system. The centripetal acceleration 

plays the role of the magnetic anisotropy, whereas the combined effect of external field 

F and the gravity field is equivalent to the external magnetic field. 
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